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ON THE USE OF SOMOFF'S THEOREM FOR THE EVALUATION OF THE 
ELLIPTIC INTEGRAL OF THE THIRD SPECIES. 

By Mr. Chas. H. Kummell, Washington, D. C. 

[CONTINUED FROM VOL. II, PAGE TJ^\ 

For the convenient computation of the elements of the next higher step, I 
propose the following algorithm. Starting with a = i , b = /3, c = y, aA<p, c cos <p, 
aAft, c cos ft, as data, compute as follows : — 

a' = %(a + c), V = \{a — e), c' = V '{ac); 

a' 4<p' = \ (aA<p + c cos <p), c' cos <p ' = i/(a' A<p' -\- b') . j/(a'A<p' — b') , 

a' Aft' = | (a J ft + c cos ft), c' cos ft' = \/{a'Aft' + V) . ^/{a'Aft' — b') , 

a' A ft" = i (aAft — c cos ft), c' cos ft' = \/{a'Ap> + b') . \/{a'Aft y — b').{i8) 

We thus obtain the new elements of the next higher step in the modular 
scale, viz. y' = c' j a', f', ft', ft\ from which (16) can be formed. Computing in 
the same manner a", b", c" from a', b', c' ; a"Af" and c" cos <p" from a'Af' and 
c' cos <p' ; a" Aft" and c" cos ft", as well as a" Aft n and c" cos ft n , from a' Aft' and 
c' cos //'; a"Aft u and <r" cos //', as well as a" Aft" and c" cos /*", from a'J// and 
c' cos ft', we obtain the elements of the equation 

n{ fy> ft y ) = n{<p>;„ t $,) - /%;-„ ^ - /%;„ #,) + /%;„ ^,). (16") 

Continuing this process until a (n) = c (n) = arithmetico-geometric mean of a 
and c, and also a^A<p^ = c (n > cos f (n >, d n) Af£ n -»' = <r (n) cos /* ( — l) ', but J^"-"' = 
cos f/"~ 1)x = o, then we have finally 

Z% y> ^) = I [(-)' //(f<" ( > ) .^" ( nl 1> ')] = * [(-)' # (tf°, M"- 1 ")]^!^) 

where s denotes the number of grave accents affixed to ft. 

The integral is now expressed in 2" terms of the same form ; viz. : — 

/7M»), M"- 1 ") = fsin/zC-^'cos 2 //^- 1 ".-^ . . ? lt if ly . , 

\n > n i J r r cos f i — sin 2 ju (n 1} sin 2 <p 

, • ,„ i« , i + sin tf w . ,i + sinH Cn_1) 'sintf> (n) , . 

= |sin//"- 1) '/— . Y M — lt — . („_!„ . y , nV (19) 

2 ' 1 — sin f^' z 1 — sin //" '' sin ^ (n) v *' 

This is the standard analytical form. To adapt it to the use of logarithmic 
tables, place 

cos / ( "- 1) ' — sin f/ n ~ 1)f sin <p (n) ; (20) 

,. 1 , 1 -f sin <p in) 1 tl _ . , , , nVv , , 

then, since <p y = ^ / ; - ^ ^ „ = ^ / tan 1 ( j + y> <•>), (21) 
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we have 

n (f{ n \ iA n - l> ') log e = sin ft* - I)( log tan £ ( J + y><»>) — log cot ^ (n ~ 1) ' 

= <r (n) sin //<" ~ :) ' if y log e — log cot £i<" " 1) '. (19') 

Computing this quantity for each limit, and combining according to (i6 (n) ), 
we obtain IJ (if u ) log e, whence FI (if u) . 

This method is not applicable to the cyclometric integral II (if , Vgi), as I 
had at first thought, although the algorithm ( 1 8) is applicable to a A (v g i)_ = a Av 
sec v, c cos (yai)_ y = c sec v in a straight line, because these and a' A (v'a/i)_ ,, 
c' cos (y\yi)_ y , and a" A (v"^,i)_ yiy , c" cos (v" 3 „i)_ y „, etc., are all real. But, 
although a' A (v s s ,i)_ y , is also real, c' cos (v'j } ,i)_ y , = i/[a /2 ^ 2 (^ / «')_ y , — £' 2 ] 
is not, because a' A (v\,i)_ , < b'. This makes the above method, which I shall 
distinguish as the first method, impracticable for the cyclometric integral, at 
least for the ascending scale. 

There is a slight defect in this method in case if is small ; for the algorithm 
( 1 8) ends with a (B> Af (n) = <r ( "> cos <p (n) , or Aw <n> = cos if <n> , which is near unity if if, 
and hence f {n) , are small. The determination of w (n) from the cosine is in that 
case not precise. 

XT • / , n/\ s i n <P' cos <p' 
Now, since sin if = ( 1 + p') a~, — 

_ a sin if' cos y>' 

, . c' sin w', c' cos #' 

we have c smip = , , . *— ; 

T a'Aip' ' 

, . . a' Aw' . , x 

. • . c' sin tc' = —. — r — ; . c smip. (22) 

r c' cos If' T x ' 

Adding this easy computation for each step of the modular scale, we have 
finally both <r (n) sin if in) and t 00 cos ip (n) from the algorithm, and f <n) can be deter- 
mined from its sine, cosine, or tangent. Moreover, since (20) requires sin if (n) 
and sin //"', this results directly. 

While this controlling computation is not essential to this method, it is to 
the second method, which is based on form (17) of Somoft's theorem. We have, 
at the second step of the scale, 

n {( f v , y) = * n w;„, #,) - 2 [f sin , %/ - v \±$2*2? ) 

j . 1 -f r sin n sin if 

— r sin a ip v + \l — — — — z - 

1 r-Ty i j — y sin// sin if 
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= 2 2 /7 (#?,„ fy) - (2c> sin p> + c sin p) <p y + I cot 2 \i! cot #, (23) 

,,„ . . c sin p . c sin w , . . 

where cos X w = 7- sin // sin f = * (24) 

,. , . , . , c' sin u'. c' sin <p' , A 

cos V = f s\n p' sin <p' = —, — -, — ; (24') 

and finally, at the limit, 

II{tp , p) = 2 n II (<p{ n \ p{ n) ) — <f y (2" - l 6 n ~ J) sin p (n ~ l) + . . . + 2c' sin p' + c sin p) 

+ 2* , - 1 /cotP"- 1) + ... + 2/cOt-^ / + /cot£A< 0) . (25) 

The first term is, of course, computed by (19')- If. however, the control (22) 
is used, the computation of the limiting step may be omitted ; for, since then 

a (n) Jyfli) _ c (n) cos ^W 

we have by (22) t (n> sin <p w = c (n ~ l) sin <p^- l \ 

C M s ; n ffio — t (— i) s i n ^c-D ; 

therefore cos A (n) = cos #" - », 

and /7 ($•>, tf») = c™ sin //»' ^ y — / cot *• A<"> 

= «:<» - >> sin //" ~ l) f y — l cot i *<• - ». 
Regarding this, we have, instead of (25), 

n (<P r t* 7 ) 
= <f (2 n - 1 c (n ~ 1) sin p (n ~v — 2"- 2 c (n - 2) sin p (n ~ 2) — ... —\2c' sin p' — c sin /i) 

— 2 ( "- 1) /cOt4A < "- 1) + 2 (n - 2) /cot^-A ( "- 2 > + . . .+ 2/cotj;' + /cot^;. 

(250 

Now, with a slight modification, this method may be used for the cyclo- 
metric integral. Placing, for the sake of uniformity, 

Py = V ( 26 ) 

then the algorithm for parameter starts with 

a dp = a d{voi)_ — a Jvsec v (27) 

and c cos p — c cos (v^z) = c sec v, (28) 

and is real in a straight line. The control (22) becomes 

c' sin 1/ = c' sin (v'n,i) , = —. — — -. . c sin p = ic' tan v', 

I \ ji' J—y' c ' COS p' r 

a' Aa' , . 

or c ' tan v = —. — - — -. . c tan v. (29) 

r cos p' 
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The formula (24) cannot be used, since X is not real. Placing, instead, 

, , c tan v .c sin . . 

tan x <0> = *-, (30) 

, , 1 + r sin a sin , , 1 + ri tan v sin tc 

-we have \l — — - — = — T - = \l —.- ; — L 

* I — y sin p. sin i 1 — ^ tan v sin y? 

_ . . 1 + ttanx (a) 
— * 1 — »tan*<°> 

= «<»>, (31) 

and therefore 

= <p y i(2 a ~ 1 tr (n - 1) tan v^- 1 ' — 2 n - V"- 2) tan »/"- 2) — . . . — 2c' tan v' — <r tan v) 

— «(2' , - 1 X < "- 1) — 2 n - 2 X (n - 2) — . . . — 2X' — X (0) ) . (32) 

There is, however, a modification of these methods required, if the modulus 
y is very small ; for, although even then the method necessarily converges, it will 
do so after a greater number of steps have been ascended, which makes the com- 
putation both tedious and less precise. If y is small then /9 is near unity. The 
required modification is therefore obtained by reducing the integral to the com- 
plementary modulus by Jacobi's imaginary transformation. We have thus 

#(? r /v) ==/7( <v ,, y ) 

***' d r 2 • 



o 



_%. J ("/?)-/? . d<!> — f tan 2 4> 

_ j i tan v cog2 v • t j^ /3) _ i3 - 1 _f tan 2 v tan 2 ^ 



Jv dtp f sin 2 (p 



/dv dip r sin' 
tanV oo?"v-^- fi .,,- 
o 1 5- sin 2 4> 



(33) 



Jv 



Since sin(j^+ Jy ,-+ v /J )_ /J = _ _ 

C08 ( J /»+ J r , '+ v /»)-/»--j9^;' 
^ ( -^ 4- -i y * + v^) _ j8 = «'r tan v ; 

we have 

/3 2 sin(J /? + J y *+ v /3 )_ /3 cos( jp-f- J r ? + v /3 )_ / jJ(J /3 + J y » + i^)..,, 



108 KUMMEIX. ON THE USE OF SOMOFF's THEOREM FOR THE EVALUATION 

= -f tan v — 5- , 

' COS*!- 

and therefore 

*fsin(j / ,+ -y' + v /J )_ j gCos(j /? + -i r * + v p )_ p ^(Jp+ -J 7 i + v^_ p 

= J I v *± i^sinV 

o ^ A #-i-^sin J (j /3 + J ) ,* + v j g)_ j8 sinV 

cos(vo — ^a)_3 
- " (<V ^ " h tan v Jv + \ ^ cos(> ; + ^_; ' * (9). 

f y tan^J/£, . , i +^tan^J/£.jtanyJy> 
= "(<VV — T' i~ + t/ i-4-tan/.J//.itan^J f 

- *<* v + *^*-* efc+ff:; (34) 

= //(^o, vg) + ^ y tan // J/* — /cot £ a», if cos <o = tan ft Jft tan 9? Jf . (35) 

The term II{<po,Va) is treated like a logarithmic integral to the modulus /9. 
We start the algorithm (18) and control (22) with the data 

a=i, c = y, b = fi; 

a J(p == a dip sec tp , £ cos <p = b sec y>, £ sin <p — - tan f ; 

a dv = a Aft sec /*, £ cos v = 3 sec //, £ sin v = -. tan ft ; (36) 

and computing the next higher step of /9, being the next lower step to 7-, the 
quantities 

a, = £ (a + £), f, = £ (« — 3), 3, = >/(*£) ; 

a, dip, = $ (a dip -{- & cos <p) = a dip, sec y>,, 

£, cos ^, = i/(a, 2 J 2 </>, — <r, 2 ) = b, sec y>,, 

... a, dip, ... 3, 

3, sin w, = -7 i-r . sin c- = -. tan to, ; 

r b, cos f , ? r 

a, dv, = \(a dv -\- b cos v) = a J ft, sec /i,, 

b, cos v, = j/(a, 2 J 2 f, — c , 2 ) = £, sec /.<,, 

£ • a ' dvi , . b, i„e.i\ 

0, sin v, = 1 . b sin v = -t tan «,. (30') 

3, cos v, 1 w ' 

Continuing in this manner until a (n) = b w = arithmetico-geometric mean of 
a and 3, we have also a (n) d<{\ n) = b {n) cos <p w or J^ (B) sec ip in) = sec y> (n) , and simi- 
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larly for ft. By (25') we have now, going back to its more fundamental form, 

= <pp (2" -1 £(,,_,) sin v {n _ l) — 2 H ~ 2 b (n _ 2) sin v (n _ 2} — ... — 2b, sin v, — b sin u) 

_ 2 »- 2/ 1 + A.-D sini, (»-i) s ' n ^(.-i) 

I — A— 1) S ' n V (» - 1) Sm <P(.« - 1) 

J. 2 . - s / - 1 + /?(»- 2) sin v (» - 2) sin ^ ( „ _ ,,) 
1 — A» - 2) sin v (n _ 2) sin ^ (n _ 2) 

+ . . . 



. 1 + 3, sin v, sin ^, , . 1 + ^ sin v sin ^ 
1 — 3, sin v, sin d>, * 1 — 8 sin v sin t£ 



= — <P Y ( 2 * " ' K - 1) tan // (n _ I, — 2" - 2 b (n _ 2) tan //<„ _ 2) — . . . — 2<*,[tan ft, —'b tan //) 

+ 2" - ' / cot $X (n _ „ — 2" " 2 / cot ^ ( „ _ 2) — . . . — 2/ cot \X, — I cot ^ (0) , (37) 
where we have placed 



, «, , £ tan u . b tan tc 

cos X m = a tan u tan tp = <- — ; 1 , 

T a .b 

in. . b, tan a, . b, tan w, 
cos ,?, = 8, tan //, tan 10, = - <-!—l 12 e tc. 



(38) 



We have also here 



*V ^ 2* (rt i - sin ft,. 






(39) 



and we recognize in this the elegant formula of Gauss for the integral of the first 
species, and <p {n) is the argument of Jacobi's & functions, usually denoted by x. 
Treating the cyclometric integral in a similar manner, we have 

n {fr v p i ) = n {<Pli i >-t l y) 

= /7(yr / 3,j / 3 + -y — /y) 

= 11(0*. - ug) + 4> tan vAv - ^ COs(v /?~^-/? 

= — II (<p B , v 3 ) — i [<p tan vAv — arc tan (tan v Av tan <p J^)] 
= — 77 (<pp, vp) — i [y y tan vAv — o), if tan o = tan v Av tan <p A<p. (40) 
For the first term the algorithm starts with the data 
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a=i, c = y> b = $; 

a A<p = a Atp sec <p, b cos <p = b sec <p , bsva. <p= -tan^>; 

aJv, bcosv, bs'mv; (41) 

and, since v is real, we can pursue the algorithm for v also sideways, by computing 
also a, Av x =^(aJv + b cos v) and following up each branch to its limit. We 
have then, according to the first method, 

/7(^,^) = ![(-)• /7([^,] 1 ,[v ( „_ ,„].)], (42) 

where 

/7/r/ n r i \ j. 1 1 / l + sin V (»-D' sin $*(«> 

/y ([^)]"[ y (-')]o=^)»"^-.)>^-^ I _ sin ^_ l ; tSin ^ it) ) 

= — (B) » sin v ( „ _ d, f y + «' arc tan (sin v ( „ _ 1)( tan <p (n) ) 
= — * (*<») sin v (n _ 1)f f y — x ( „ _,),), 
if tanz ( „_ 1)f = sinv ( „_ I) , tan^ w . (43) 

By the second method we have 

= ^ /3 (2"- : /5 { „_ 1) sin v ( „_ 1} — 2"- 2 £ ( „_ 2) sinv ( „_ 2) — . . . — 2b, sin v, — b sin v) 

_ 2 » -2 / I +A»-l) S ' nt '(n-l) Sin l(,,-l) 

1 — A» - 1) sin v <» - 1) sin ^(» - 1) 

, 2 «-3{ l + /?(n-2)Sinv ( „_ 2) SU1 ^ ( „ _ 2) 
1 — /9 ( „ _ 2) sin v ( „ _ a) sin </> (n _ 2) 

+ • • • 

. 1 + & sin v, sin <p, , , 1 + /? sin v sin ^ 
1 — /?, sin v, sin ^, * 1 — /} sin v sin ^ 

= — f r «'(2"-M ( „_ 1) sin !/(„_„ — 2"- 2 3 (n _ 2) sinv ( „_ 2) — . . . — 2<*,sin v, — bsm v) 

+ t(2"- x X in _ l) — 2"- 2 X ( „_ 2) — . .. — 2X, — X (0) ), (44) 

where tan * (s) = /9 W sin v (s) tan f w = *c> »in v (t) A> tan y w (4S) 

I have thus shown that the algorithm which forms the main part of the 
computation is applicable in all cases of the integral of the third species which 
need to be considered. The practical use of the methods will be shown in a 
future article. 



